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We present a set of dynamical equations based on Ashtekar's extension of the Einstein equation. 

The system forces the space-time to evolve to the manifold that satisfies the constraint equations or 

the reality conditions or both as the attractor against perturbative errors. This is an application of 

the idea by Brodbeck, Frittelli, Hiibner and Reula who constructed an asymptotically stable (i.e., 

O^ ' constrained) system for the Einstein equation, adding dissipative forces in the extended space. The 

^\ , obtained systems may be useful for future numerical studies using Ashtekar's variables. 

,-H ■ PACS numbers: 04.20.Ex, 04.20.Fy, 04.20.Ha 
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"*^ \ I. INTRODUCTION 

m ; 

^^ . Ashtekar's extended formulation of general relativity [|l[ has many advantages. By using his special pair of variables, 

_. 1 ' the constraint equations which appear in the formulation become low-order polynomials, and the theory has the correct 
►^ , form for gauge theoretical interpretation. These features suggest the possibility for developing a nonperturbative 
fvq ' quantum description of gravity. When we apply his formulation to classical dynamics, however, we have to impose 
\^ , the reality condition additionally to the system, since Ashtekar's variables will not produce a real-valued metric in 
^^ ■ general. 

>0 ' Fortunately, it was shown that the secondary condition of the reality of the metric will be automatically preserved 

during the evolution, if the initial data satisfies both primary and secondary metric reality conditions . If we impose 
the reality condition on the triad (triad reality condition), then we have additional conditions that can be controlled 
by a part of a gauge variable, triad lapse ^q (defined later) [^. Therefore the reality conditions are controllable, and 
^ ■ we think that applying the Ashtekar formulation to dynamics is quite attractive, and broadens our possibilities to 
I \ attack dynamical issues. 

^ ' It is, however, also the case that preliminary numerical simulations of the spacetime using Ashtekar's variables 

ZV\ show that the system will not normally recover real-valued spacetime if we relax the metric reality condition locally 
J> . during the evolution B. Therefore, we desire a system that is robust for controlling both the reality conditions and 
the constraint equations for stable long-term integration. 

In this article, we propose new dynamical systems based on Ashtekar's variables, which satisfy this desire. That 
C^ ' is, even if numerical data gives us a truncated solution which violates the constraints or the reality conditions during 
steps of time evolution, the system forces the spacetime to evolve to a manifold that satisfies the constraint equations 
and the reality conditions against these small deviations. This is an extension of the idea by Brodbeck, Frittelli, 
Hiibner and Reula (BFHR) g] who constructed an asymptotically stable system (i.e., it approaches to the constraint 
surface) for the Einstein equation. (A similar effort can be found also in B.) BFHR introduced additional dynamical 
variables. As, which obey dissipative dynamical equations and which evolve the spacetime to the constraint surface of 
general relativity as the attractor in the extended spacetime. 

Recently, we have obtained a set of equations of motions for Ashtekar's canonical variables (i?* , A'^) in a symmetric 
hyperbolic form W (sec also |H,p[). We here present first that a set of constraint equations also forms a symmetric 
hyperbolic system in its evolution equation. This is the Ashtekar version of the work by Frittelli pC| , and is already 
shown by Iriondo et al. 0. Our version, however, is based on the fixed inner product throughout evolution. We then 
show that a dissipative system for the combination of {E]^,A1) and constraints also forms a symmetric hyperbolic 
system, following the procedure of X-system by BFHR. We expect that this system has an attractor in the constraint 
surfaces. 
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We next extend this system also to one that has an attractor in the real-valued surfaces. Since the dynamical 
system that we have obtained in a symmetric hyperbolic form requires the triad reality condition ||7|, our purpose is 
to construct a system which asymptotically evolves into the triad-real-valued manifold. We show this is available by 
applying the same technique. 

In this article, we discuss only the case of a vacuum spacetime, but including matter is straightforward. 

II. ASHTEKAR FORMULATION 

We start by giving a brief review of the Ashtekar formulation and the way of handling reality conditions. We also 
describe a symmetric hyperbolic system that was obtained in M . 

A. Variables and Equations 

The key feature of Ashtekar's formulation of general relativity Q is the introduction of a self-dual connection as 
one of the basic dynamical variables. Let us write the metric 5^^ using the tetrad Ej^ as g^^, ~ Ej^E'lrju [|. Define 
its inverse, i?j, by E'^ := E^g^^^rju and we impose /?° = as the gauge condition. We define S0(3,C) self-dual 
and anti self-dual connections ^^° :— (jJ^^^ ^ («/2)e°bc w^°, where w^'^ is a spin connection 1-form (Ricci connection), 
Lol^'' := E^'^'Vfj^E^ . Ashtekar's plan is to use only the self-dual part of the connection +^" and to use its spatial part 
'^Ai as a dynamical variable. Hereafter, we simply denote +^" as „4". 

The lapse function, N, and shift vector, A^', both of which we treat as real-valued functions, are expressed as -Eq = 
{1/N, —N'^/N). This allows us to think of Eq as a normal vector field to S spanned by the condition t = x^ =const., 
which plays the same role as that of Arnowitt-Deser-Misner (ADM) formulation. Ashtekar treated the set (i?* , A'^) 
as basic dynamical variables, where E'* is an inverse of the densitized triad defined by _E* := eE"*, where e := detEf 
is a density. This pair forms the canonical set. 

In the case of pure gravitational spacetime, the Hilbert action takes the form 

S ^ j d'x[{dtA'i)El + {^/2)NKEiF^,^'''c 

-e'AN-N'Ft,Ei+A^„V,K], (2.1) 

where N := e^^N, F° :— 29[^^°i —ie°'bc Al^ii^l is the curvature 2-form, A is the cosmological constant, ViEi :— diEl — 

'^lE^, where e„fe := eabcEfE]. 



ieab'AlEi, and e^ = detE^ = (detEf)^ is defined to be detE^ = (l/6)e"''^e,jfcE^£;^E;^, where e^Tc := €abcEfE''E^ 



^ijk • — ^ ^ijk ri- 

Varying the action with respect to the non-dynamical variables N, iV* and Aq yields the constraint equations, 

Ch ■■= (tl-lYK KHf^j - a det E « 0, (2.2) 

Ca« := -F^Ei « 0, (2.3) 

Caa ■■= V,El « 0. (2.4) 
The equations of motion for the dynamical variables (E* and Af) are 

dtK = -^T^ii^^'aNEiEl) + 2V,{N^=EJ) 

+iA''oeab' El, (2.5) 

dtA^ = -le'^^c NElF^^ + N^F^^ + V.A^ + ANE^, (2.6) 



^We use ij.,u — 0, ■ ■ ■ ,3 and i, j = 1, • • ■ , 3 as spacetime indices, while I,J = (0), ■ ■ ■ , (3) and a,b = (1), ■ ■ ■ , (3) are SO{l, 3), 
SO{3) indices respectively. We raise and lower fj,, v, ■■■ hy g'^" and gfj,i, (the Lorentzian metric); /, J, • ■ ■ by i]^'^ — diag( — 1, 1, 1, 1) 
and riu; i,j, • ■ ■ by 7'-' and 7;^ (the 3-metric); a,b,- ■ ■ by S'^'' and Sab- We also use volume forms eabc- ^abc^'^'"^ = 3!. 
When {i,j,k) = (1,2,3), we have e^^fc = e, e^jk = 1, e"-'* = e~\ and e*^* = 1. 
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where P^-Xf := djXl' - leab^A^Xl^ for X^^ + XI' = 0. 

B. Reality conditions 

In order to construct the metric from the variables {El^,Af,N,N'^), we first prepare the tetrad Ej as Eq = 

{l/eN, ~N'/eN) and E^^ = (0, El/e). Using them, we obtain the metric gi^" such that g^''' := E'^E'^jr/". 

This metric, in general, is not real-valued in the Ashtekar formulation. To ensure that the metric is real-valued, we 
need to impose real lapse and shift vectors together with two metric reality conditions; 

ImiKEn = 0, (2.7) 

W'^ := Reie'''"'E^El'VkE^^^) = 0, (2.8) 

where the latter comes from the secondary condition of reality of the metric Im{9f (-E*i?^°)} = [^, and we assume 
deti^>0 (see @). 

For later convenience, we also prepare stronger reality conditions, triad reality conditions. The primary and sec- 
ondary conditions are written respectively as 

K ■■= MK) = 0, (2.9) 

and Im(9ti^^) = 0. (2.10) 



Using the equations of motion of E^, , the gauge constraint J|2.4| ), the metric reality conditions (2.7), ( |2.8| ) and the 
primary condition (^), we see that ( ^.1C| ) is equivalent to ^ 

Re(^g) = d^{N)E''' + {l/2e)E^NE^''djEl 

+N'Re{Af), (2.11) 

or with un-densitized variables, 

Re(^g) = d.,{N)E"' + N' Re{A^). (2.12) 

From this expression we see t hat t he secondary triad reality condition restricts the three components of the "triad 



lapse" vector Aq. Therefore ( 2.11 ) is not a restriction on the dynamical variables (£'* and Af) but on the slicing, 
which we should impose on each hypersurface. 

Throughout the discussion in this article, we assume that the initial d ata o f (El ^A° ) for evolution are solved so as 



to satisfy all three constraint equations and the metric reality condition (2/7) and ( |2.§| ). Practically, this is obtained, 



for example, by solving ADM constraints and by transforming a set of initial data to Ashtekar's notation. 

C. A symmetric hyperbolic form 

We say that the system is a first-order (quasi-linear) partial differential equation system, if a certain set of (complex) 
variables Ua (a = 1, • ■ • , n) forms 

dtUc,^M'^^c.{u)diup+Mc,{u), (2.13) 

where M. (the characteristic matrix) and J\f are functi ons o f u but do not include any derivatives of u. If the 



characteristic matrix is a Hermitian matrix, then we say (2.13) is a symmetric hyperbolic system 



For a pair of Mq = {El,A^), a symmetric hyperbolic system is obtained by modifying the right-hand-side of the 
dynamical equations using appropriate combinations of the constraint equations. The final form of the system is 
written as 

M{E, E)^"^'^ = ie'^'^'NEiY' + N^fH"^, (2.14) 

M{E,Ay''^'^ = MiA^Ey^^'^ = 0, (2.15) 

MiA.A)^"^'^ = iNie^^-'El-l^' - e'^^-'El-^i' 

_g-2^»ag&cd^j^i^ - e'^e'^'^^EiiEiE'^ 

"•^E'^EIE^^) + N^S^^-i'^ (2.16) 



_g-2gacd 



where A4{*, *) means a block component of the characteristic matrix as 



At, 



MiAEY''^, M{A,AY\^' 



A) 



(2.17) 



where = means that we have extracted only the terms which appear in the principal part of the system. The inner 
product of a set of the variables is 



{{K^AmK.At)) = l^J5'''KEi + f'^abAlA]- 



(2.18) 



We note that this symmetric hyperbolic system is obtained under the assumption of the triad reality condition, 
together with gauge conditions, A^ = A^N^ and diN — 0. 



III. ASYMPTOTICALLY CONSTRAINED SYSTEM 



Frittelli |l0| showed that the propagation of the constraint equations in the standard ADM system of the Einstein 
equation forms a symmetric hyperbolic system. This fact suggests that a small violation of the constraint equations 
such as a truncation error in numerical simulation does n ot b eh ave in a fatal way immediately. 

Similarly, we can show the set of constraint equations (2.2), ( ^.3[ ) and (2.4), forms a symmetric hyperbolic system 
in its evolution equations. The principal part of the time derivatives of Ch, Cmi '■— ^Cmi and Cca become 




^ Q 



iNea'^'E'^ + N^Sa" J 




(3.1) 



which forms a Hermitian matrix under the inner product rule of 

{{CH,CMi,CGa)\{CH,CMi,CGa)) '■ — 



(3.2) 



(C) _ 



We note that non-principal parts of the dynamical equations of a set of Uq = {Ch,Cm,Cg) include the terms of 
Ua and Ua ■ These facts suggest that all constraint equations have a well-posed feature. Iriondo et al [^ present a 
similar result, but our definition of the inner product does not include any coefficients. We also remark that all other 
occurences of the inner products throughout this article also do not include any coefficients (i.e., obey the normal 
index notation). Thus we omit to express the inner product hereafter. 

Following the BFHR procedure Iq], we next construct a dynamical system which evolves the spacetime to the 
constrained surface, Ch ~ Cmi ~ Cca ~ as the attractor. We introduce new variables (A, Ai, Aa), as they obey the 
dissipative evolution equations 



dfX = ai Ch ~ Pi A, 
dtXi = a2 Cmi — 02 Xi, 

dtXa = a^CGa - hXa, 

where a^ ^ (allowed to be complex numbers) and j3i > (real numbers) are constants. 

If we take Ua — {El, Af, A, Ai, Aa) as a set of dynamical variables, then the principal part of (|3 
written as 

dtX ^ -iaie'>^''El&MA]), 
dtX, ^ a^i-eSlEl + e5iEi]{diA''j), 
dtXa = a^diEl 

The characteristic matrix of the system u„ does not form a Hermitian matrix. However, if we modify the right- 
hand-side of the evolution equation of {El^,At), then the set becomes a symmetric hyperbolic system. This is done 



(3.3) 
(3.4) 
(3.5) 

can be 

(3.6) 

(3.7) 
(3.8) 



by adding asY^idiXa) to the equation of dtEl, and by adding iaie" c"^ E^ &^{di\) + a2(-e7'™i?f + e(S™i?'")(9/A™) to 
the equation of dtAI- The final principal part, then, is written as 



/5\ 

At 
A 

VaJ 



/MiE^EYa 



I bi 



\ asSa'S' 



b:l 







ib 



'iaieb'"^E^E 




a2e(5ri^^ - ^^i^D 





he^/E^El^ a2e{S^E^'' - 7'"^"; 











/ 



A'' ' 
A 

\xj 



(3.9) 



Clearly, the solution (E'^,^", A, Ai, Aa) = {El, Af, 0,0,0) represents the original solution of the Ashtekar system. If 
the As decay to zero after the evolution, then the solution also describes the original solution of the Ashtekar system 

in that stage. Since the dynamical system of Mq , (3.9), constitutes a symmetric hyperbolic form, the solutions to 
the A-system are unique. BFHR showed analytically that such a decay of As can be seen for As sufficiently close to 
zero with a choice of appropriate combination of as and /3s, and that statement can be also applied to our system. 
Therefore, the dynamical system, (p^), is useful for stabilizing numerical simulations from the point that it recovers 
the constraint surface automatically. 



IV. ASYMPTOTICALLY REAL- VALUED SYSTEM 



We next extend the system, ( |3.9| ), to the one that also has an attractor in the real- valued surfaces. Since the 
dynamical system of Ua requires the triad reality condition to form a symmetric hyperbolic system, our purpose is 
to construct a system which asymptotically evolves into the triad-real-valued manifold. 

In order to obtain such a system, we add two more new variables, A*^ and A^, which satisfy the evolution equations 






(4.1) 
(4.2) 



corresponding to the secondary metric reality conditions, ( [2.8|) , and the primary triad reality conditions, (2.9), re- 
spectively. The equation (4.1) is necessary to complete this system since the term W^^ appears in the non- principal 
part of the equation dtUl. Thus, these two equations will guarantee the consistency with the secondary condition. 
The principal terms of (4.1) and (4.2) become 



dtX''9^aiRe{e'^'''E';Ei^dkEl^), 
dtXl - 0, 



(4.3) 
(4.4) 



under the assumption of the triad reality condition. 

A set of variables Ua — {El,Af,X,Xi,Xa,X^^,Xl^), then, forms a symmetric hyperbolic system if we further 

modify the equations for dfEl similarly to the case of (pM. The final set of equations in the matrix form can be 
written as 



Af 
X 

X^ 
Xa 



( 



aiRe{e^^'^E[E^^5'l) 
\ 



X(3.9) 



a,Re{ea^dE[EfJl.) 
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f Er\ 
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X 
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Xm 

Xb 
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0/ 




VAr/ 



(4.5) 



where A^([3.9|) denotes the matrix in the right-hand-side of (3.9). 

Clearly, the solution {El,At, X, Xi, Xa, A'^ , A^) = {E^, Af, X, Xi, Xa, 0, 0) represents the original solution of (3^). The 
same discussion in the previous section can be applied also here. Therefore, we expect that the system (4.5) controls 
the violation of the triad reality condition during the time integration. 



We remark that a reduced set of the variables Ua 



(DR) 



{El, Af, X'^^ ,Xa) does not work for the purpose of controlling 



the reality condition, since the secondary condition of the reality requires the constraints to be satisfied. 



V. CONCLUDING REMARKS 

We showed a set of dynamical equations, which has the constraint surface as its attractor, by introducing new 
additional variables that obey dissipative equations of motion. Based on BFHR's analytical proof ^ , we expect that 
this set of equations is robust against a perturbative error of the constraint equations. Thus, the system may be 
useful for future numerical studies with its stability property. 

We also showed an advanced set of equations that has its attractor also in the real-valued surface. Since our 
symmetric hyperbolic system of the original Ashtekar's variables requires the reality condition on the triad, the new 
system is designed as such a way. The same above discussion can be applied to this advanced set, and we expect the 
asymptotically real-valued feature in its evolution. 

The problem of these systems might be that they require many additional variables. From a view point of numerical 
applications, this claim would not be so serious a problem, as we see a success of a dissipative maximal slicing condition 
|ll| ('K-driver' in the literature). Actually, the A-system of the Einstein equations was already tested and confirmed 
to work appropriately in numerical applications at least in one dimensional space evolution models |l2| . Therefore we 
expect our system also shows the desired asymptotic behaviours. We are in preparation of presenting such a numerical 
result. We are also trying to reduce the number of the variables in order to find out clear geometrical meanings of 
our A-system. We will report these efforts elsewhere. 
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